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ABSTRACT 

We present observations at 3.6, 4.5, 5.8, 8, 24, 70 & 160^m of NGC 205, the dwarf elliptical 
companion of M31, obtained with the Spitzer Space Telescope. The point sources subtracted images 
at 8 and 24^m display a complex and fragmented infrared emission coming from both very small 
dust particles and larger grains. The extended dust emission is spatially concentrated in three main 
emission regions, seen at all wavelengths from 8 to 160/im. These regions lie approximately along 
NGC 205's semi-major axis and range from ^100 to 300 pc in size. Based on our mid-to-far infrared 
flux density measurements alone, we derive a total dust mass estimate of the order of 3.2 x IO^Mq, 
mainly at a temperature of ^ 20K. The gas mass associated with this component matches the predicted 
mass returned by the dying stars from the last burst of star formation in NGC 205 (~0.5 Gyr ago). 
Analysis of the Spitzer data combined with previous 1.1mm observations over a small central region 
or "Core" (18" diameter), suggest the presence of very cold (T~ 12K) dust and a dust mass about 
sixteen times higher than is estimated from the Spitzer measurements alone. Assuming a gas to dust 
mass ratio of 100, these two datasets, i.e. with and without the millimeter observations, suggest a 
total gas mass range of 3.2 x 10^ to 5 x IO'^Mq. 

Subject headings: galaxies: dwarf — galaxies: individual (NGC 205) — galaxies: ISM — infrared: 
galaxies — infrared: ISM — ISM: dust, extinction — ISM: structure 



1. INTRODUCTION 

The galaxies of the Local Group, because of their prox- 
imity and our ability to carry out high spatial resolution 
studies in nearby systems, provide us with powerful tem- 
plates to understand the chemical, morphological and 
kinematical characteristics, plus star formation history, 
of more distant galactic systems. The Local dwarf galax- 
ies (LDGs), in particular, have been subject to numerous 
recent studies of their stellar populations and star forma- 
tion (for summaries sec for example, Hodge 1989, Mateo 
1998, Grebel 2005). Dwarf ellipticals follow a different 
"fundamental plane" than true ellipticals (Ferguson & 
Binggeli 1994), nevertheless they do share many common 
characteristics, in particular the presence of a mostly old 
stellar population (Mateo 1998). One of these local dwarf 
systems is NGC 205, a low surface brightness dwarf ellip- 
tical companion of M31. NGC 205 is interesting because 
of its very conspicuous dark clouds, which were detected 
in some of the earliest photographic plates of M31 and 
its surroundings (Baade 1944; Hodge 1973) (NGC 205 hes 
~ 36' NW from the M31 nucleus). These visually dark 
clouds reveal the presence of an interstellar medium in a 
stellar population dominated by old stars, that is, in an 
elliptical galaxy. 

There are several indicators for an intermediate and 
young stellar population, as well as for gas enrichment, 
in NGC 205. A handful of blue and UV-bright stars 
have been detected in its nucleus (Baade 1951; Hodge 
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1973; Bertola et al. 1995), with ages ranging from 
~ 5 x 10* yr to a few 10^ yr. In terms of metallicity, 
the location of its giant branch implies a mean value of 
[Fe/H] = -0.85±0.2 (Mould, Kristian & Da Costa 1984), 
higher than for most LDGs, which have values ranging 
from -1.1 to -2.2 (Mateo 1998). Moreover, a stellar 
population study of the resolved bright near-IR stars 
in NGC 205 (Davidge 2003) indicates that the bright- 
est asymptotic giant branch stars formed a few tenths 
of a Gyr ago. All these indicators suggest an episode of 
star formation in NGC 205, which took place a few 10* 
years ago, and a recent replenishment of the interstellar 
material by this population of evolved stars. Although 
there is no consensus on the process which triggered such 
a burst, gravitational tidal interaction of NGC 205 with 
M31 has always been an appealing mechanism because 
it simultaneously explains the twisted surface brightness 
isophotes (Hodge 1973; Choi et al. 2002), the tidal trail 
of blue metal poor RGB stars (McConnachic et al. 2004; 
Ibata et al. 2001) and the kinematically distinct behavior 
and morphology of its interstellar medium, particularly 
of the HI gas (Young & Lo 1997). 

The ISM in NGC 205 has some peculiarities. Based 
on the last burst of star formation ~ 5 x 10* yr ago 
(Wilcots et al. 1990; Bertola et al. 1995), one can esti- 
mate the amount gas injected back into the ISM, which 
is - 2.6 X IO^Mq (following Faber & Gallagher 1976). 
This mass is about three times larger than that mea- 
sured by Sage, Welch & Mitchell (1998) of 7.2 x lO^M© 
(from CO (H2) plus HI, scaled by a factor 1.4 to include 
Helium). This difference in masses led to the idea that 
gas is being removed by the young stars in NGC 205, 
either by stellar winds or supernova explosions (Welch, 
Sage & Mitchell 1998). Recently Haas (1998) carried out 
photometric observations with the Infrared Space Obser- 
vatory (ISO) between 120 and 200/im that, combined 
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with a millimctric measurement of the central 18" (Fich 
& Hodge 1991), yielded large amounts of cold dust with 
mass 0.2 - 1.2 x lO^'M© at T ~ 10 K. Assuming a gas- 
to-dust mass ratio of 100, the cold dust traces nearly 
10^ Mq of gas at center. There is vast observational 
evidence for "warm" and "cool" dust in elliptical galax- 
ies based on mid/far infrared data (see e.g. Jura et al. 
1987; Knapp et al. 1989; Ferrari et al. 2002; Temi et al. 
2004; Xilouris et al. 2004). The initial IRAS measure- 
ments were limited to detection of temperatures of ^ 30 
K traceable by the lOOfiui fluxes. This led to the pos- 
sibility of missing a colder component with nearly 90% 
of the dust mass at lower temperatures (T ~ 15 — 20K) 
which is radiating at longer wavelengths (> 100/im) (see 
e.g. Devereux & Young 1990). The longer wavelength 
far infrared observations (> 100/im) combined with mil- 
limeter and sub-millimeter measurements have confirmed 
this picture in several early-type galaxies, where the bulk 
of the dust mass is cold (~ 20 K) and nearly an order of 
magnitude larger than previously estimated (Wiklind & 
Henkel 1995; Temi et al. 2004). 

This study presents new mid/far infrared images of 
NGC 205 obtained with the Spitzer Space Telescope 
(Spitzer). The higher sensitivity and spatial resolution 
enable us for the first time to look directly into the 
detailed structure and properties of the dust clouds in 
NGC 205 and to evaluate their dust properties and con- 
tent. The outline of the paper is as follows. In Section 2, 
we describe the Spitzer observations and the data reduc- 
tion. The Spitzer miages are presented in Section 3, along 
with a description of the mid/far-IR emission morphol- 
ogy and a comparison with other tracers of the ISM. Ad- 
ditional image processing is discussed in Section 4. Dust 
mass measurements and gas mass estimates are derived 
in Section 5 using two different methods. In Section 6, 
we discuss the results with some of the previous mass es- 
timates and theoretical expectations. We conclude with 
a summary of our findings in Section 7. 

2. OBSERVATIONS 

The mid-IR observations were obtained with IRAC, 
the infrared camera on board Spitzer (Fazio et al. 2004) 
on 20 January 2005 as part of the GOl program to 
map M31 (PI: Barmby, Program ID: 3126). The region 
around NGC 205 was observed with four dithered 30 sec 
frames per sky position; the IRAC field of view (FOV) 
is 5'. The images obtained with channels 1 and 2 (at 
3.6 and 4.5/im) are used mostly to determine the stellar 
contribution to the overall spectral energy distribution of 
NGC 205. The channels 3 and 4 (at 5.8 and 8/im) contain 
the most information about the dust properties; chan- 
nel 4 in particular, with a passband centered at 7.9/im 
includes some of the strongest Polycyclic Aromatic Hy- 
drocarbons (PAH) emission features known (e.g. bands 
at 7.7 and 8.6/im); these are considered some of the best 
tracers of very small particles (see e.g. Li & Draine 2001, 
Draine & Li 2001 and references therein). The Basic Cal- 
ibrated Data produced by the SSC pipeline (Version 11) 
were combined with the SSC MOPEX software to pro- 
duce final mosaics with pixel size 1.2" (all channels) and 
spatial resolutions of less than 2" FWHM. 

The far-IR data were obtained with the multiwave- 
length photometer MIPS onboard Spitzer on 25 August 
2004 (Program ID 99). These observations were taken 
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Fig. 1.— IRAC 3.6, 4.5, 5.8 & 8/im images of NGC 205. The 
FOV of each image is 15' X 15' with N up and E to the left. 



fiis figure was to big for astropin, it is JPG now 



Fig. 2.— IRAC Sfira MIPS 24, 70 & 160Atm images of NGC 205. 
The FOV of each image is 5' X 5' with N up and E to the left. 



as part of the large M31 mosaic consisting of 7 scan legs, 
between 0.75 and 1.25 degree long, with a 148" over- 
lap between scan legs. The scanning was done at the 
medium scan rate covering a total area of 1° x 3° ap- 
proximately oriented along the major axis of M31, with 
a total observing time of 17.1 hrs. The NGC 205 obser- 
vations themselves have a depth of 84 sec per pixel. The 
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Fig. 3. — A three color image of NGC 205 as seen by IRAC 
at 3.6 (blue), 5.8 (green) and 8/xni (red). The image emphasizes 
the dust clouds distribution seen at longer wavelengths over a 15' 
square region, with N up and E to the left. 



Fig. 4. — The three main emission regions — "North" , "Center", 
"South"- as well as the "Core" and "Total" regions (as defined in 
Table 1) are shown ovcrlayed on the IRAC Sfim imago of NGC 205. 
The FOV is 5' x 5' with N up and E to the left. 



details of the observations and reductions are given by 
Gordon et al. (2006). In the three MIPS bands a back- 
ground subtraction has been performed, that removes 
some time dependent artifacts along the scan legs. In 
the 70 and 160/^m images some stripping remains due to 
a non-linear time response of the Ge:Ga arrays, however 
their overall effect on the flux density around NGC 205 
is very small (~ 1 — 2%), and this is included in our un- 
certainties (see Section 4.1). The final 24, 70, and 160 
micron images have spatial resolutions of 6, 18 and 40" 
FWHM, respectively (Ricke et al. 2004). The mosaic 
pixel scales are 2.49" at 24//m, 9.85" at 70/^m, and 16.0" 
at 160/um. 

The Spitzer data arc complemented with ISOCAM ob- 
servations of NGC 205 from the ISO archive. The archive 
data (Off Line Processing VIO.O) was re-reduced using 
the CIA (Cam Interactive Analysis V4.0) software to 
improve mainly on the deglitching, following the rest of 
the standard reduction steps (Coulais & Abergel 2000; 
Blommaert et al. 2003). The data set includes images at 
14.3^m covering a 204" x 216" field of view and with 6" 
per pixel. All the measurements were performed with 
background subtraction, and were compared with the 
global estimates of Xilouris et al. (2004) for consistency. 

For both IRAC and MIPS wc used a 15' x 15' section 
of the M31 maps centered on the position of NGC 205 
[a(2000) == 00''40"22.1^ 5(2000) = 41°41'07.1"]. This 
region is large enough to encompass most of visual size of 
the galaxy of 9' x 19' (25 B mag arcsec'^) (Nilson 1973), 
and certainly all the mid/far IR emission (Rice et al. 
1988). We adopt a distance to NGC 205 of 815 ± 35 kpc 
(Mateo 1998; McConnachie et al. 2004) throughout this 
work, i.e. at that distance, spatial structures 1" across 
measure ~ 3.9 pc. This tentatively puts NGC 205 behind 
M31, which is at a distance of 780 kpc (Rich et al. 2005; 
McConnachie et al. 2005). 

3. MORPHOLOGY OF INFRARED EMISSION 



Figure 1 and 2 show the images of NGC 205 taken 
with IRAC at 3.6 & 4.5/^m (stellar emission dominant, 
see Table 4) and at 5.8, 8, 24, 70 & 160/im (dust emission 
dominant, sec Table 4). Figure 3 displays a three color 
image of NGC 205 where the dust clouds at 5.8 and Sfim 
stand out within the stellar background detected by the 
shorter wavelength bands. The higher angular resolution 
images at 8 and 24/im display a spatially complex and 
fragmented extended emission. We expect that the main 
contribution to the 8/im emission arises from aromatic 
molecules and at 24/im from submicron-sizcd grains (Li 
& Draine 2002; Xilouris et al. 2004). The overlay of 
these two wavelengths (Fig. 4) shows the close spatial 
relationship between these two dust distributions. How- 
ever, the 8/im image shows more extended emission than 
at 24/im (see Fig. 5). 

The morphology of the diffuse emission is complex and 
clumpy; nevertheless we can clearly identify three spa- 
tially distinct emission regions: the largest region to the 
south with a size of ~ 76" x 37" (~ 300 pc x 146 pc), 
the central region with a size of ~ 49" x 28" (~ 193 pc x 
110 pc), and the smallest region to the north with a size 
~ 33" X 31" (~ 130 pc x 123 pc) (see Fig. 4 and Table 1). 
Even at the lower resolution of the 70 and 160/im bands 
(see Fig. 2), where the bulk of the emission is produced 
by dust emission from the larger grains (~ 0.1-0.5/im, 
Desert, Boulanger & Puget 1990; Li & Draine 2001), one 
can easily identify the same three large regions seen at 
shorter wavelengths. We have included also a "Core" 
region with a diameter of 18. "4 for which 1.1mm obser- 
vations are available (Fich & Hodge 1991). 

The 14.3/im ISOCAM observations, however, cover a 
smaller FOV and the final mosaic has some vignetting 
that leaves an effective FOV of ~ 192" radius. Although 
our mosaic is a bit larger than that presented by Xilouris 
et al. 2004 (Fig. A.l), the 'South" and "Total" regions 
are 3% and 6% smaller at 14.3/J,m, respectively, than at 
all other wavelengths (see Table 3). Therefore, the emis- 
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TABLE 1 
NGC 205 Emission Regions 

Region^' RA (J2000) Dec {J2000) Area 

h m s o / // /2 

North 40 24.00 41 41 49.9 0.24 

Center 40 22.29 41 41 09.4 0.30 

South 40 23.75 41 40 10.3 0.55 

Core 40 22.29 41 41 09.4 0.07 

Total 40 22.38 41 40 52.3 4.36 
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^elliptical aperture as shown in Fig. 4 



"This figure vv/as to big for astroph, it is JPG novv/ 



Fig. 6. — The HI column density derived by Young & Lo (1997) 
are shown overlayed on the IRAC Sfim grayscale image of NGC 205 
where point sources have been removed. The HI density contours 
range from 2 xlO^^ cm~^ to 3.8 xlO^" cm~^ in intervals of 4 
xW^^ cm-2. The FOV is 5' X 5' with N up and E to the left. 



Fig. 5. — IRAC 8/xm grayscale image of NGC 205 where point 
sources have been removed to reveal more clearly the extended dif- 
fuse emission. The image is displayed at high contrast to enhance 
the faint emission surrounding the bright central region. The FOV 
is 15' X 15' with N up and E to the left. 



sion not included directly in these two regions, given the 
very low surface brightness at 14.3/im, is small. 

A comparison of the mid/far IR data with other trac- 
ers of the ISM, such as HI (Young & Lo 1997), shows 
the same morphology (sec Fig. 6), indicating that dust 
and gas are very well mixed in NGC 205. Based on the 
HI velocity data, it is clear that the HI gas rotates in 
the North-South direction with a velocity gradient of 
42 km s"^ over 900 pc (Young & Lo 1997). However, 
there is no clear one-to-one correspondence between the 
HI column density peaks and velocity, and therefore, the 
three distinct IR emission regions are not necessarily dy- 
namical entities. The observed dust morphology is simi- 
lar at all wavelengths from 8 to 160/im (within the limits 
of varying angular resolution) and also similar to the HI 
morphology seen in the 21-cm line. The data is consis- 
tent with an ISM in which the atomic gas and dust are 
well mixed and therefore we see no evidence for signifi- 
cant stratification of the ISM in NGC 205. 

4. ADDITIONAL DATA PROCESSING 

A reliable estimate of the IR and diffuse emission in 
NGC 205 requires an estimate of the foreground point 
source contamination arising either from the Milky Way 
(MW) or M31. Point sources were extracted in the cen- 
tral 5' X 5' region of NGC 205 in the 5.8, 8.0 and 24/xm 



images using the source extraction software StarFinder 
(Diolaiti et al. 2000). Spitzer [8]-[24] and [5.8]-[8] col- 
ors were derived for these point sources and compared 
to the colors of a sample of Galactic AGB and K Gi- 
ant stars templates in the mid-infrared observed with 
the Short Wavelength Spectrometer (SWS) on board of 
ISO. Using ISA? (ISO Spectral Analysis Package Ver- 
sion 2.0)^ it is possible to determine the flux density 
within each of the Spitzer's imagcr/photometer bandpass 
over the SWS spectra which cover a wavelength range of 
2.4-45.2^m. The point sources in the central region of 
NGC 205 (Fig. 7 top, open squares) have colors simi- 
lar to those of AGB stars (filled triangles). However, 
K Giant stars in the field are clearly identifiable due to 
their much "bluer" [8]- [24] color and, consequently, we 
conclude that none are present in the central region of 
NGC 205. 

The foreground contamination at 24/im, from M31 
point sources, where dusty AGBs contribute most of their 
flux, is estimated using six control fields. Each field has 
dimension 15' x 15' to match that of NGC 205 and is 
located a comparable distance away from the center of 
M31 (Table 2). To determine the number of AGB stars in 
NGC 205, we count these stars in the central 15' square 
region and estimate the number of foreground stars by 
averaging their numbers in the six control fields. Based 
on the results shown at the bottom of Figure 7, we esti- 
mate that about half of the point sources in the NGC 205 
15' X 15' field are from foreground contamination due ei- 
ther to the MW or M31. 

We estimate that in the smaller "Total" region cen- 
tered on NGC 205 (as defined in Table 1), about a third 
of the point sources are from foreground contamination. 
Given a mean flux density of 0.6 mJy for contaminating 
point sources (sec Fig. 7, bottom), we estimate a flux 
density contamination of ~ 2 mJy. This number is be- 

^ ISAP is available at|http://www.ipac.caltech.edu/iso/isap/isap.html| 
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Fig. 7. — Top: Color-color diagram of a sample of K Giant stars 
(filled diamonds) , AGB stellar templates {filled triangles), and the 
point sources extracted in the central 5' of the NGC 205 field (open 
squares). Bottom: Histogram showing the number of point sources 
as a function of 24/xm flux in the central 5' of the NGC 205 field 
(solid line), for the average of the six control fields (see Tabic 2; 
dashed line), and the remaining distribution for NGC 205 after the 
average from the control fields has been removed (dotted line). 



low our flux density measurement errors at 24/xm (see 
Table 3) and therefore no further processing, i.e. point 
source removal, is done to the 24//m image. 



TABLE 2 
NGC 205 Pointing and Control Fields 


Field 


RA (J2000) 

h m s 


Dec (J2000) 

o / // 


NGC 205 
control field a 
control field b 
control field c 
control field d 
control field e 
control field f 


40 22.10 
36 38.65 
44 54.11 
43 36.11 
46 24.97 
47 38.73 
42 54.21 


41 41 07.1 
40 37 47.0 
40 51 03.5 

40 31 49.4 

41 10 13.8 

41 28 52.5 

42 19 51.6 



4.1. Aperture Correction of Extended Sources 

We measure the IRAC 3.6, 4.5, 5.8, & Sfim and the 
MIPS 24, 70, & leO^m fluxes in the spatially resolved 
regions common at all wavelengths, as depicted in Fig- 
ure 4 and listed in Table 1. The IRAC and MIPS data are 
combined with previous observations taken with ISO at 
U.SniJi (Xilouris et al. 2004) and JCMT at 1.1mm (Fich 
& Hodge 1991). The flux uncertainty for each region in- 
cludes the combined effect of the systematic errors due 
to the post processing of the images, plus that of the flux 
density calibration of the instruments (Fazio et al. 2004; 
Rieke et al. 2004). The flux uncertainty, therefore, was 
estimated by measuring the residual flux density away 
from the galaxy in the background subtracted images 
and adding an absolute flux calibration uncertainty of 
10% for the IRAC and MIPS 24/im measurements, and 
of 20% for the ISO U.3^lm. and MIPS 70 & leO^im fluxes. 
The SfiiR fluxes are corrected (multiplied) by the effec- 
tive aperture correction factor of 0.944 (3.6^m), 0.937 
{A.dfim), 0.772 (S.S^im), and 0.737 (8//m) (Reach et al. 
2005). The corrected flux densities and uncertainties are 
given in Table 3. 

TABLE 3 
Photometry of Emission Regions in NGC 205 



Wavelength 


North 


Center South 


Core 


Total 


fim 




Flux Density in mJy 




1.22 


430±2 


130±7 110±6 


46.2±2 


940±47 


1.65 


550±3 


160±8 40±7 


56.8±3 


116±58 


2.16 


380±2 


110±6 100±5 


39.1±2 


820±41 


3.6 


21±3 


59±6 54±7 


23.5±2 


438±52 


4.5 


13±2 


35±4 33±5 


14.9±1 


266±33 


5.8 


15±3 


33±5 39±9 


14.2±1 


279±62 


8 


12±1 


26±3 26±3 


14.5±1 


148±15 


14.3 


7±2 


15±3 25±5 


5.9±1 


101±23 


24 


10±4 


18±5 26±10 


5.8±1 


141±23 


70 


104±21 


181±36 301±60 


61.7±12 


1428±290 


160 


473±279 


713±363 1490±692 


253±45 


8761±4720 



4.2. Stellar Photospheric Component of Diffuse 
Emission 

The stellar photospheric component of the diffuse emis- 
sion coming from NGC 205 itself must be removed be- 
fore the analysis of the dust component. We conserva- 
tively assume that, in addition to the 2MASS J, H, and K 
bands, all the light in the 3.6 and 4.5 IRAC bands is pho- 
tospheric in origin. The J, H, K, IRAC 3.6 and 4.5 filter 
response functions were multiplied with a PEGASE stel- 
lar model (Fioc & Rocca-Volmerange 1997) with an age 
of 350 Myr (see Fig. 9) and integrated over wavelength 
yielding fluxes that can be compared with the observed 
fluxes. The age was picked to roughly coincide with the 
(uncertain) epoch of recent star formation in NGC 205 
(see eg. Butler & Martincz-Delgado 2005; Davidge 2005). 
There is no diagnostic power in the age assumed; we used 
synthetic spectra from Simple Stellar Population (SSP) 
models only to estimate the stellar contribution to the 
Spitzer bands, as simple power laws cannot reproduce the 
J, H, K colors. Scalings were derived for each band indi- 
vidually and averaged together to get the final scaling for 
the photospheric component. In all cases, the J through 
4.5/im fluxes are consistent with being completely stellar 
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in origin. To estimate the photosphcric contribution to 
the flux density at longer wavelengths, we multiplied and 
integrated over wavelength the relevant IRAC and MIPS 
bandpasses with the PEGASE model and subtracted the 
result from the observed flux density to yield an estimate 
for the total dust flux density. For the relative contribu- 
tions of the diffuse photospheric emission to the observed 
fluxes, see Table 4. The age we use for the PEGASE 
models (within a factor of 3) does not significantly alter 
the dust fits; indeed, PEGASE models between 0.1—1 
Gyr result in very similar estimates for the stellar con- 
tamination in the Spitzcr bands, especially for A > Sfim. 

TABLE 4 
Stellar Photospheric and Dust Flux Density Contributions 

Region Wavelength TotaP Dust Stellar'' Dust Fraction 
/im nijy mjy 



North 



3.6 

4.5 

5.8 

8.0 

14.3 

24.0 

70.0 

160.0 



20.8 
12.6 
15.2 
11.7 
6.9 
10.1 
104.2 
472.9 



0.0 
0.0 
6.4 
6.6 
5.2 
9.5 
104.2 
472.9 



1.0 
0.60 
0.42 
0.25 
0.08 
0.03 
0.00 
0.00 



0.00 
0.00 
0.42 
0.56 
0.75 
0.94 
1.00 
1.00 



Center 3.6 


58.6 


0.0 


1.0 


0.00 


4.5 


35.4 


0.0 


0.60 


0.00 


5.8 


32.6 


7.2 


0.43 


0.22 


8.0 


25.5 


10.7 


0.25 


0.42 


14.3 


14.9 


9.8 


0.08 


0.66 


24.0 


17.6 


15.7 


0.03 


0.89 


70.0 


181.5 


181.5 


0.00 


1.00 


160.0 


712.9 


712.9 


0.00 


1.00 


South 3.6 


54.1 


0.0 


1.0 


0.00 


4.5 


32.8 


0.0 


0.60 


0.00 


5.8 


38.8 


15.5 


0.40 


0.40 


8.0 


25.7 


12.1 


0.23 


0.47 


14.3 


25.5 


20.9 


0.09 


0.82 


24.0 


25.8 


24.0 


0.03 


0.93 
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''The measured flux density (photospheric plus dust contribu- 
tions) 
''Stellar flux density normalized with respect to that at 3.6/im 



5. MASS OF DUST AND GAS IN NGC 205 

With these new Spitzer observations, it is possible to 
derive the mass of dust and estimate the gas mass in 
NGC 205 using different methods and compare the re- 



sults, for the total emission as well as within the three 
main emission regions. 

In the first method, we use a single-temperature fit to 
the 70 and 160/im data to determine the dust tempera- 
ture, assuming a modified black body emissivity function 
of the form: 



F{X) ^ \-^- 



2hc^X- 



,hc/kXT _ I 



(1) 



We adopt an emissivity coefficient (3 = 2 (valid be- 
yond 70/.tm; cf. Draine 2003). Assuming a gas-to-dust 
mass ratio of ~100, representative of the solar neighbor- 
hood and the inner disk of the Milky Way (Hildebrand 
1983; Devereux & Young 1990), the dust temperature 
and emissivity can be used to estimate the mass of dust 
and gas. This method has been extensively used to de- 
termine dust masses in galaxies, but it has a few draw- 
backs (Devereux & Young 1989; Dale & Helou 2002) and 
should be taken as an approximation. A fit to the spec- 
tral energy distribution (SED) has been shown to provide 
better estimates to the dust masses for each dust grain 
component (see below). 

A reliable ratio of the far-IR fluxes is a necessary ingre- 
dient to get the best temperature estimate. The angular 
resolution of the 70/.tm image was reduced to match that 
at 160/xm by convolving the image with the 160/^m ker- 
nel, derived from the PSF model generated using the 
TinyTim software (Krist 2002). The IRAF routines ge- 
omap and geotran were then used to match the pixel 
scale and astrometry, preserving the flux in the trans- 
formation. Each pixel is assigned a temperature based 
on the best fit to the modified blackbody function. The 
resulting temperature map is shown in Figure 8. The 
temperature is very uniform, ranging from ~ 14 — 17 
K in the innermost 5' that encompasses the three main 
emission regions. In the central rectangular region of size 
2.1' X 2.9' (or an area of 6.3'^, slightly larger than the 
"Total" region as defined in Table 1), we calculate the 
dust mass from the 160/im emission using the following 
equation: 



Md = S,DyKdB,,T 



(2) 



where D is the distance to the galaxy of 815 kpc and Kd 
is the mass absorption coefficient for which we assume a 
value of 6.73 m^ kg^^ at 70/im and 1.2 m^ kg^^ at 160/im 
(Li & Draine 2001). The total dust mass in this region, 
derived from the sum of the individual mass estimate for 
each pixel, is 6.1 ± 1.2 x W^Mq, yielding a gas mass of 
6.1 X IO^Mq for a gas-to-dust mass ratio of 100. 

Alternatively, we can determine dust masses by fitting 
the total dust emission SED. The SED was measured 
between 5.8 and 160/(m, using IRAC, ISO 14.3/(m and 
MIPS data. We measure the fluxes in the spatially re- 
solved regions common to all wavelengths, as depicted in 
Figure 4 and listed in Table 3. 

The dust emission in NGC 205 was modeled by an 
equation of the form 



FdustW — y ^CiKi{X)B\{TD^i) 



(3) 



where the the sum extends over i dust components. 
Ci = Md.il D"^ 1 D as in Eq. 2, and Ki is the mass absorp- 
tion coefficient of the i*'* dust component. For NGC 205, 
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Fig. 8. — Dust temperature of NGC 205 within the central 
innermost 5' with 24^m eontours shown for levels 0.3, 0.35, 0.4, 
0.45 and 0.5 MJy/sr. The temperatures are based on the 70 and 
160/xm flux ratio. The 24 and 70^m images have been convolved 
to the 160/xm resolution. 



the relatively large 5.8, 8, and 14.3^m fluxes require a hot 
component which we ascribe to PAH molecules. While 
the 70 and 160/im emission could be fit with a single 
temperature component, including the 24^m and 1.1mm 
fluxes requires additional temperature components. To 
satisfactorily fit the full SED of NGC 205, we adopt a 3 
component model: warm and cold silicates (a ~ 0.1/im) 
and PAHs, the former to reproduce the 24-1100/im emis- 
sion and the latter, the 5.8-24^m emission. Mass absorp- 
tion coefficients for astronomical silicates were computed 
from Mie theory using the dielectric functions of Laor & 
Draine (1993). Cross sections for the PAH molecules 
were taken from Li & Draine (2001). As the canonical 
PAH spectrum (e. g. NGC 7027, Werner et al. 2004) 
exhibits no features beyond 20^m, we re-computed the 
PAH cross-sections (and mass absorption coefficients) 
leaving off the last three terms of Li & Draine's Eq. 11. 
The PAH component is included here for completeness of 
the fit, but we do not use the temperature or mass esti- 
mates as the PAH emission is a stochastic process rather 
than an equilibrium process as assumed in Eq. 3. 

With only 6 measurements and 6 parameters 
(C'i,2,3,Tdi 2 3)1 ^ formal fit is not fully constrained and 
we proceed by making a detailed search of parameter 
space. With our assumptions about the components 
present in the observed SED, we searched a limited set of 
parameters for dust mass (through Ci) and temperature 
for each component. We assumed temperature ranges of 
200-550 K, 25-150 K, and 8-25 K for the PAH, warm, and 
cold silicates, respectively. At each set of temperature 
points, a similar grid of mass scalings {Ci) was explored 
for each component. After an initial, coarse grid search, 
a finer temperature and mass scaling grid was defined for 
each component and the procedure repeated. We iterate 
this procedure (typically only two refinements of the grid 
were needed) until a good fit in the x^ sense was found. 



We estimated uncertainties in the derived quantities by 
examining the distribution of x^ about the minimum for 
all input parameter sets (typically 1 x 10^ — 1 x 10^ points 
in the six dimensional parameter space). While the re- 
sults of this procedure are not unique, it yields represen- 
tative estimates of both the dust mass and temperature, 
given the input assumptions. In all cases the dust mass 
is dominated by the cold component and should be con- 
sidered as a weak lower limit. Owing to the lack of data 
beyond 160/im, we cannot formally rule out much larger 
dust masses. However, several pieces of evidence argue 
against dust masses at the upper end of the formal error 
bars. First, where we do have sub-mm data ("Core" re- 
gion), the dust mass is consistent with much lower upper 
"limits" on the mass and second, the upper error limits 
would imply very large fluxes at sub-mm wavelengths, 
which is not what is observed (Fich & Hodge 1991). 

TABLE 5 
Mass and Temperature Estimates from SED Fits 



Cloud 



North 



Warm Silicates 



1-2: 
53 



iiMo 



2.2 



Cold Silicates 



X 10^ Mq 



K 



17.5 



+5 



K 



Center 



2.1+2Af0 

,+11 

-6 



53"' 



K 



2.9_i 7 X IO^Mq 
17+1 K 



South 



5.3 
51 



+ 18 
3 



K 



Core'' 



Core*^ 



53t| K 

1.9ti Mq 

48tl K 



1.2J:o g X 10"' M0 
15.81:^ K 



l-llag X IO^Mq 



17.i 



:;+3 



11.6+S K 



Total 



6.1 

58 



TT3- 

— 5 
+25 



Afo 



3.2 



-1.8 

18.1 



W^Mp 



+2 



''SED without the 1.1mm flux density measurement 
''SED includes the 1.1mm flux density measurement 

The best SED fits to the data for each region, which 
produce the dust masses and temperatures given in Ta- 
ble 5, are shown in Figure 9 and 10. Using the same 
gas-to-dust mass ratio of 100, we infer a gas mass for the 
"Total" region of 3.2 x lO^Af0. This value is nearly one 
half of the 6.1 x IO^Mq mass derived from the 70 to 160 
fiux ratio method. Summing up the total dust mass in 
the three distinct large regions using the best fit values, 
we estimate a gas mass of 1.7 x lO^M©, i.e. smaller than 
that of the "Total" region. This result suggests that a 
larger reservoir of gas/dust with a comparable gas mass 
surrounds the three regions, as first suggested by the HI 
(Young & Lo 1997) and CO (Welch, Sage & Mitcheh 
1998) observations. 

The "Core" region has a detection at 1.1mm (Fich & 
Hodge 1991), and adding this value to its SED and com- 
puting a new fit, results in a dust mass of 1.8 x W^Mq, 
i.e. approximately a factor of sixteen larger than the mass 
estimated without the 1.1mm point (1.1 x 10'^ Mq), and 
at a colder temperature, 11.6 K (compared to 17.8 K). 
One can speculate that if this colder dust component has 
a similar distribution across our selected regions, preserv- 
ing the same mass gas-to-dust ratio, then one predicts a 
gas mass several times larger, ~ 5 x lO^il/©. This mass 
determination is consistent with the fact that adding sub- 
millimeter measurements to far-IR observations of early- 
type galaxies usually yields larger cold dust mass esti- 



10. OF 




FIR Mapping of NGC205 

10.0 F 



10 100 



1000 




10 100 

X [/u.m] 



1000 



10. OF 




10 100 

X [/xm] 



1000 



Fig. 9a. — Infrared SEDs of the "North" [top left) and "Center" 
{top right) region of NGC 205 (see Table 1). The diagrams show the 
SED fit to the IRAC, li.Sfim ISO, and MIPS measurements (Table 
3). The solid points are the data (stellar photospheric component 
removed) and the open points are the model multiplied with the 
passband and integrated over wavelength. The dotted line is the 
PAHs component, the dashed line is the warm silicates component, 
the dot-dashed line is the cold silicates component and the solid 
line is the total of all components. 



mates (as much as a factor ten) than the far-IR mea- 
surements alone (Devereux & Young 1990; Tcmi ct al. 
2004). 

6. DISCUSSION 

One of the open issues in understanding the ISM of 
NGC 205 has been its "missing ISM mass", i.e. the 
apparent lack of gas in comparison with the expected 
mass return to the interstellar medium by evolved stars 
(Welch, Sage & Mitchell 1998). By determining the 
dust mass and selecting a reasonable gas-to-dust mass 
ratio (e.g. 100), one can estimate the amount of gas in a 
straightforward manner. 

The dust mass estimate for the "Core" of NGC 205 
from Fich & Hodge (1991) using IRAS data and including 
1.1mm observations is Mdust = 1 — 3 x 10"^ Mq, depend- 
ing on the adopted value for the dust emissivity. The re- 
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Fig. 9b. — Infrared SEDs of the "South" {bottom left) and 
"Total" {bottom right) region of NGC 205 (see Table 1). 



cent mass dust determination in the same region by Haas 
(1998), adding ISO data between 120 and 200Atm to the 
SED, is Mdust = 0.88 x IO^Mq without the 1.1mm data 
and Mdust = 1-18 x W^Mq including it (for a emissivity 
exponent /3 = 2), i.e. a factor of four higher than Fich 
& Hodge's upper limit. Using the Spitzer data alone, we 
estimate for the "Core" region Mdust = 1.1 x 10^ Mq. If 
we include the 1.1mm data, the mass estimate increases 
to Mdust = 1.8 x lO^Af©. This latter value is similar to 
that obtained by Haas, except that our dust temperature 
is 11.6 K rather than 8 K. Therefore the mass of gas, for 
a gas-to-dust ratio of 100, for NGC 205 "Core" region is 



Mga. 



gas 



1.8 X lO'^M, 



0- 



Using only our FIR Spitzer data for the SED, we de- 
termine a "Total" amount of dust at 18 K of Mdust = 
3.2 X IO^'^Mq, which implies a mass of gas, Mgas = 3.2 x 
lO^Af0. This later value is four times larger than ISM 
mass determination of Mjsm = 7.2 x IO^Mq by Sage, 
Welch & Mitchell (1998), and in perfect agreement with 
the theoretical gas mass estimate oi Mgas = 2.6 x 1O^M0, 
generated by the injection of gas by evolved stars over 
0.5 Gyr. This comparison suggests that the mass deter- 
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Fig. 10.— Infrared SEDs of the "Core" region of NGC 205, as 
defined in Table 1. The diagram shows the SED fit to the IRAC, 
14.3/.im ISO, and MIPS measurements (Table 3) plus the 1.1mm 
data point measured by Fieh & Hodge (1991). The symbols and 
lines are defined in Fig. 9. 



mination by Sage, Welch & Mitchell (1998) is indeed a 
lower limit to ISM mass, and that essentially over our 
"Total" region one can account for all the ISM expected 
in NGC 205 from the highly evolved stars since the last 
starburst. 

If the mass of dust in NGC 205 is indeed at a few 
times larger than previous estimates, then we believe 
this provides a more coherent picture of its recent star 
formation history. The global evolutionary path of the 
ISM in NGC 205 was summarized by Welch, Sage & 
Mitchell (1998). They pointed out that any reasonable 
efficiency during the last starburst in NGC 205 implies 
a more massive ISM. Their idea is straightforward: if 
the burst transformed ^^ 1.4 x IO^Mq of gas into stars 
(the Wilcots et al. 1990 value adjusted to include low 
mass stars < IMq) with a 10% efficiency, then NGC 205 
must have had at least ^ 1.2 x IO'^Mq of gas. Our dust 
mass estimates could easily account for at least ~ 1/3 
of such a gas reservoir, and nearly all if the gas temper- 
ature were ^ 14K. Finally, a burst of ~ 1.4 x 10^ Mq 
would yield ~ 90 stars in the brightest magnitude of the 
AGB according to the fuel burning theorem of Renzini 
& Buzzoni (1986) [assuming a mass to light ratio of 0.2, 
the same as that of NGC 1866, a 10* year old star clus- 
ter (Fischer et al. 1992)]. This is of the same order of 
magnitude as we see in Figure 7 (bottom). 

7. SUMMARY 

New IRAC and MIPS observations of NGC 205 give 
a better assessment of the total dust mass in this pro- 
totypical LDG system. The Spitzer data have been 



complemented with ISOCAM archival data to create a 
well defined spectral energy distribution in the 3.6 to 
160/im wavelength range. The higher sensitivity and an- 
gular resolution of IRAC and MIPS provide a detailed 
view of the morphology of the interstellar medium in 
NGC 205. At 8 and 24/im the spatial structure of the 
warm dust is clearly resolved. Although the 8 and 24/xm 
images resemble the drawings of the dust clouds inferred 
by Hodge (1973) using isodensity measurements of the 
ground based visual photographic plates (c.f. Welch, Sage 
& MitcheU 1998, Fig. 1; Young 2000, Fig. 1), these new 
IR data display a more complex structure. To facilitate 
our analysis of the color temperature and dust masses, we 
identified three main regions. As shown in Figure 2, the 
dust clouds within these regions are bright and conspic- 
uous at the longer wavelength too, having sizes ranging 
from ~ 100 to 300 pc, and dust masses of ~ IO^-IO^Mq. 
A larger "Total" area encompasses the three regions, and 
covers most of the bright 160/^m dust emission. 

The derived dust masses, based on the 70 to 160^m 
flux ratio and a best fit to the SED, imply gas masses 
(assuming a standard gas-to-dust mass ratio of 100 (cf 
Young 2000)), that are consistent with each other and 
larger (3.2 — 6.1 x IO^Mq) than previously estimated from 
CO and HI observations. This indicates that a significant 
contribution to the dust and gas masses is coming from 
a fainter component surrounding the brightest emission 
regions. Adding the 1.1mm detection to our measure- 
ments of the "Core" region of NGC 205, increases by a 
factor of sixteen the estimate of the cold dust in that 
region. Our gas mass estimate in NGC 205 is a factor of 
four larger than previously detected, but it could be still 
larger if the distribution of the very cold dust (~ 12 — 14 
K) spreads out to the three selected regions. If this is 
the case, the gas mass would be sixteen times the "To- 
tal" (Table 5), i.e. 5 x lO^M© for a gas-to-dust mass ratio 
of 100. Overall the gas estimates are therefore consistent 
with the predicted mass return from dying stars, based 

yr ago. 
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